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Abstract 


This  investigation  was  undertaken  as  a  part  of  a  larger  effort  to  determine  the 
implications  of  using  coherent  turbulent  structures  or  turbules  to  calculate 
acoustic  scattering  from  turbulent  media,  a  process  that  has  been  suggested 
as  a  possible  non-line-of-sight  means  of  detecting  enemy  assets  on  the 
battlefield.  I  present  an  analytical  solution  to  a  correlate  of  the  incompress¬ 
ible  Navier-Stokes  equation.  The  correlate  equation  is  for  the  enstrophy,  the 
square  of  the  vorticity.  The  solution  is  an  expression  for  the  time  history  of 
the  enstrophy  for  a  particular  choice  of  the  initial  velocity  distribution. 
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1.  Introduction 

This  investigation  was  undertaken  as  part  of  a  larger  effort  to  determine 
the  implications  of  using  coherent  turbulent  structures  to  calculate  acoustic 
scattering  from  turbulent  media,  a  process  that  has  been  suggested  as  a 
possible  non-line-of-sight  means  for  detecting  enemy  assets  on  the  battle¬ 
field.  I  have  published  a  number  of  papers  on  the  scattering  of  acoustic  sig¬ 
nals  from  atmospheric  turbulence  on  the  battlefield,  as  well  as  a  new  way 
of  modeling  turbulence  [1].  My  new  modeling  method  is  called  the 
Turbule  Ensemble  Model  (TEM).  The  TEM  method  consists  of  representing 
the  turbulence  region  as  a  collection  of  turbules  (eddies)  of  prescribed  sizes 
and  random  locations.  A  turbule  is  an  isolated  inhomogeneity  in  either  the 
temperature  or  the  velocity.  Methods  for  determining  the  parameter  values 
required  to  complete  the  calculation  of  the  signal  scattered  into  a  shadow 
zone  have  been  reported  elsewhere  [1]. 

As  methods  of  determining  the  parameter  values  have  improved,  reexami¬ 
nation  of  assumptions  in  light  of  theory /experiment  comparisons  has  pro¬ 
ceeded  apace.  One  of  these  assumptions  is  that  Taylor's  frozen  turbulence 
hypothesis  holds.  Previous  work  [1]  indicates  that  an  interplay  of  the  scat¬ 
tering  geometry  and  the  location  of  the  shadow  zone  boundary  determines 
the  turbule  sizes  that  must  be  considered  in  any  particular  scenario.  For 
example,  a  high  boundary  reduces  the  contribution  of  large  turbules 
because  the  scattering  pattern  of  large  turbules  is  concentrated  in  the  for¬ 
ward  direction.  This  means  that  the  scattered  signal  passes  above  the 
detector,  which  is  usually  located  near  the  ground.  The  persistence  of 
smaller  turbules  is  shorter  than  that  of  larger  turbules,  as  can  be  seen  from 
dimensional  analysis  of  the  relevant  flow  parameters  [2].  As  is  characteris¬ 
tic  of  dimensional  analysis  results,  the  constant  multiplier  to  be  used  with 
the  rate  so  determined  is  unknown.  Since  the  rate  constant  is  usually 
included  in  an  exponent,  the  effect  of  not  knowing  the  multiplier  more  pre¬ 
cisely  can  give  improper  indications  as  to  the  source  of  any  discrepancy 
between  predicted  and  measured  results.  As  the  scenario  model  is  honed 
ever  nearer  to  the  correct  representation  of  the  experiment,  improved  accu¬ 
racy  in  parameter  estimation  becomes  more  important. 

The  original  purpose  of  the  work  reported  here  was  to  come  up  with  a  bet¬ 
ter  estimate  of  the  time  interval  over  which  Taylor's  frozen  turbulence 
hypothesis  may  be  relied  upon.  As  detailed  in  section  6,  this  objective  was 
not  reached,  in  a  practical  sense,  because  the  analysis  is  for  an  isolated 
velocity  distribution.  The  mathematical  development  of  the  investigation 
is  reported  here  because  it  may  be  useful  in  research  on  decay  times  and 
perturbation  effects  as  the  Reynolds  number  transitions  from  the  laminar 
to  the  turbulent  regime.  A  short  synopsis  of  the  mathematical  development 
follows. 

Since  the  interest  at  this  time  is  on  the  evolution  of  a  velocity  turbule,  the 
Navier-Stokes  (N-S)  equation  is  the  logical  place  to  begin  such  an  investi¬ 
gation,  because  the  equation  relates  the  rate  of  change  of  fluid  velocity  to 
other  flow  quantities.  The  convective  acceleration  portion  of  the  N-S 
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equation  (the  velocity  equation)  can  be  transformed  by  a  vector  identity 
into  the  cross  product  of  the  velocity  and  the  vorticity,  plus  the  gradient  of 
the  square  of  the  velocity.  When  the  curl  of  the  resulting  equation  is 
performed  (assuming  constant  density  and  viscosity),  all  gradient  terms 
drop  out,  leaving  a  second  equation  (the  vorticity  equation),  which 
comprises  two  terms  in  the  vorticity  and  one  other  term.  The  latter  term  is 
the  curl  of  the  velocity-vorticity  cross  product.  For  a  particular  initial 
velocity  distribution,  namely  the  cross  product  of  the  angular  velocity  and 
the  position  vector  multiplied  by  a  Gaussian  envelop  function,  the  latter 
term  turns  out  to  be  parallel  to  the  velocity.  Since  the  vorticity  associated 
with  this  velocity  distribution  is  perpendicular  to  the  velocity,  the  dot 
product  of  the  vorticity  with  the  latter  term  is  initially  zero.  Thus,  a  third 
equation  (the  enstrophy  equation)  formed  by  the  dot  product  of  the  vortic¬ 
ity  and  the  vorticity  equation  has  the  interesting  property  that  the  cross- 
product  term  is  initially  zero.  The  term  enstrophy  is  defined  elsewhere  [3] 
as  the  time  average  of  the  squared  vorticity.  Enstrophy  as  used  here  refers 
to  the  square  of  the  vorticity.  In  a  further  step  in  the  analysis,  a  trial  time- 
dependent  enstrophy  distribution  expression  was  derived  by  solving  the 
enstrophy  equation,  assuming  the  cross  product  term  is  zero.  A  surprising 
result  is  that  the  time-dependent  vorticity  and  velocity  distributions 
inferred  from  the  time-dependent  enstrophy  distribution  produced  a 
cross-product  term  in  the  enstrophy  equation  that  was  identically  zero  for 
all  time.  These  velocity /vorticity  expressions  are  not  solutions  to  the  veloc¬ 
ity  equation  or  to  the  vorticity  equation.  The  expression  for  the  time  his¬ 
tory  of  the  enstrophy  is  unusual  to  the  extent  that  the  standard  definition 
of  the  time  constant  is  no  longer  appropriate. 

The  organization  of  the  remainder  of  the  report  is  as  follows:  In  section  2, 
transformation  of  the  N-S  velocity  equation  into  an  equation  for  the  vortic¬ 
ity  is  outlined.  Section  3  contains  an  analysis  of  the  characteristics  of  the 
N-S  convective  derivative  term  (as  transformed  into  vorticity)  for  the  par¬ 
ticular  initial  velocity  distribution  that  I  have  chosen.  I  show  that  the  dot 
product  of  the  convective  derivative  term  and  the  vorticity  is  zero  and, 
therefore,  the  dot  product  of  the  vorticity  with  the  vorticity  equation  elimi¬ 
nates  this  term  from  immediate  consideration.  Such  a  dot  product  without 
the  convective  derivative  term  produces  a  differential  equation  in  the 
enstrophy  only.  In  section  4,  Fourier  transform  (FT)  and  Laplace  transform 
(LT)  theories  are  used  to  solve  for  the  time-dependent  enstrophy  field.  In 
section  5,  time-dependent  velocity  and  vorticity  fields  are  inferred  from 
the  enstrophy  field  and  inserted  into  the  original  enstrophy  equation.  This 
step  confirms  that  the  inferred  time-dependent  fields,  like  the  initial  fields, 
produce  a  zero  convective  derivative  term,  and  that  time-derivative  and 
Laplacian  terms  counterpoise  each  other  so  that  the  inferred  fields  are  a 
consistent  solution  to  the  enstrophy  equation.  Section  6  contains  calculated 
results  and  recommendations  for  future  work.  Section  7  is  a  brief  sum¬ 
mary.  Appendix  A  is  a  short  review  of  applicable  FT  theory  and  notation. 
Appendix  B  contains  details  of  the  equation-solution  procedure.  Appendix 
C  contains  details  of  the  confirmation  process.  Appendix  D  contains  proof 
that  the  time-dependent  velocity  field  inferred  in  section  5  is  not  a  solution 
to  the  N-S  equation. 
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2.  The  Navier-Stokes  Equation  for  Vorticity 


The  N-S  equation  is  given  [4]  as  the  following: 


du 

dr 


■f1] 

Vp  + 

1 

V2u  + 

IpJ 

{ p) 

\P  3PJ 

V(V-u) 


(1) 


Following  the  reference,  this  equation  for  the  velocity  u  and  the  pressure  p 
may  be  converted  into  an  expression  in  which  the  gradient  terms  are 
absent.  Using  the  following  standard  vector  formulas  [5]  given  in  the  first 
line  of  equation  (2),  the  second  term  of  equation  (1)  may  be  transformed  as 


(l/2)Vw2  =  u  x  (V  x  m)  +  (w  -  V)u 

du 
dr 


PJ 


Vp  - (1/2)V u2  +  mx (V x u)  + 


(V 

Ip 


,  c 

v2«  +  —  +  — 

P  3p 


(2) 


V(V  •  u)  . 


If  the  density  p  and  the  viscosity  coefficients  r\  and  f  are  constant,  taking 
the  curl  of  both  sides  results  in  the  following  relation: 


V  x  (VII)  =  0;  if  II  is  a  scalar  function 
d  =  V  x  u 

(3) 

—  =  Vx(«xd)  +  vV2d- 
dr 

The  first  expression  in  equation  (3)  is  a  standard  vector  formula  [5].  The 
second  expression  defines  the  vorticity  d.  The  kinematic  viscosity  coeffi¬ 
cient  is  v  =  p/p,  where  tj  is  the  dynamic  viscosity  coefficient.  In  equation 
(3)  the  functions  (u,  d)  are  functions  of  the  space  coordinates  yj  ,(j  =  l,  2, 3), 
whose  units  are  typically  the  meter,  and  the  time  r,  whose  unit  is  typically 
the  second.  As  noted  in  the  introduction,  the  interest  of  this  report  is  in  the 
dynamics  of  an  atmospheric  velocity  turbule.  A  velocity  turbule  has  a  char¬ 
acteristic  linear  dimension  a  and  a  characteristic  angular  velocity  magni¬ 
tude  Qq.  These  parameters  are  conveniently  used  to  convert  equation  (3)  to 
the  unitless  coordinates  x,  the  unitless  time  f,  the  unitless  velocity  v,  and 
the  unitless  vorticity  c,  according  to  the  relations  in  the  following  equation: 


—  =  vxx(5xc)  +  KV2c  , 
dt 

where  x,  =  y,/«;  t  =  ao^;  v  =  u/(O.0a);  c  =  d/Q0; 

d  d  (4) 

V*  =  flV  ;  K  =  v/(fl2Qo);  dt  =  T-;  di  =  —■ 
y  dt  dxi 

The  notation  used  above  is  summarized  as  i  =  1,  2,  3  =  subscript  index 
(repeated  subscripts  imply  summation),  x\  =  space  variables,  x,  =  space 
unit  vectors,  and  x  =  space  vector  =  XjXj. 
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The  parameter  K  is  the  reciprocal  of  the  typical  Reynolds  number.  The 
symbol  K  is  used  rather  than  the  symbol  Re~l  for  writing  economy  The 
subscript  x  will  be  dropped  from  the  V  operator  in  subsequent  develop¬ 
ments  and,  where  appropriate,  the  shorthand  notation  for  partial  deriva¬ 
tives  indicated  in  the  last  line  of  equation  (4)  will  be  used  for  the  same 
reason. 

3.  The  Nature  of  the  N-S  Equation's  Second  Term 

The  second  term  referred  to  in  the  N-S  equation  is  the  curl  of  the  cross 
product  of  v  and  c  in  equation  (4)  for  which  the  symbol  s(x,t)  is  reserved. 
The  class  of  solutions  sought  is  for  rotationally  symmetric  velocity  distri¬ 
butions.  A  member  of  this  class  is  the  distribution  defined  in  the  equation 

u(y)  =  [Qxy]exp[-(y;/fl)2]  .  (5) 

In  equation  (5),  Q  is  an  angular  velocity  vector  and  a  is  the  characteristic 
turbule  size.  No  loss  of  generality  is  suffered  if  u(x)  is  further  particular¬ 
ized  by  choosing  Q  to  be  aligned  with  the  y$  axis  as  follows: 

^  =  V  3^0 

u0(y)  =n0[-y1y2  +y2y1]exp[-(yi/a)2]. 

Converting  to  unitless  quantities  using  equation  (4)  results  in  the  following 
velocity  distribution: 

x  =  y/a ;  v0  =  u0/(aQ.0) 

(7) 

v0(x)  =  (-x^2  +  x2Xj)exp[-(x,)2] . 

The  corresponding  vorticity  function  follows  from  equations  (3)  and  (4): 

c0(x)  =  2[x1x1x3  +  x2x2x3  +x3(1  ~x\  -x2)]expHx,)2]  •  (8) 

The  dot  product  of  the  velocity  and  vorticity  is  seen  to  be  zero,  so  the  two 
are  perpendicular. 

The  crossproduct  of  the  velocity  and  the  vorticity  is 

y0xc0  =  2[xixi(l  -xi  -  xl)  +x2x2(l  -xi  -  xi)  -  x3(xi  +x2)x3]exp[-2(x,)2]  .  (9) 

Take  the  curl  of  equation  (9)  to  obtain  the  second  term,  which  is 

So  =  4:[xiX2  -  x2xi]x3exp[-2(x,)2]  •  (10) 

The  second  term  so  is  seen  to  be  oppositely  directed  to  Vq.  This  means  that 
the  dot  product  of  the  vorticity  and  So  is  zero.  The  equation  to  be  solved  is, 
therefore,  the  following: 

■^3/C2  =  c-Vx(yxc)  +fCc-V2c  .  (11) 

This  equation  is  referred  to  as  the  transport  equation  for  instantaneous 
enstrophy  [6]  and  was  used  for  a  different  purpose  there.  When  (vq,  cq)  are 
substituted  in  equation  (4),  the  second  term  is  eliminated.  The  procedure 
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to  be  followed,  then,  is  to  solve  equation  (11)  without  the  v  cross  c  term  for 
the  initial  vorticity  distribution  of  equation  (8),  and  then  substitute  the 
result  in  equation  (11),  including  the  v  cross  c  term,  to  confirm  that  the 
answer  leads  to  a  zero  second  term  for  all  time,  and  that  the  other  terms 
combine  to  make  equation  (11)  true  for  all  time. 


4.  Transformation  of  the  N-S  Equation 

FT  theory  will  be  used  in  section  4  to  solve  equation  (11),  less  the  v  cross  c 
term,  acknowledging  in  the  process  that  this  term  is  the  very  essence  of 
turbulence.  The  applicable  theory  is  outlined  in  appendix  A,  along  with 
the  notation  employed,  and  the  definition  of  the  FT  of  a  function. 

The  convention  used  for  field  functions  is:  lower-cased  letters  are  in  the 
space  domain  and  upper-cased  letters  are  in  the  FT  or  wavenumber 
domain.  Additionally,  (v,  V)  =  velocity  field  function  (vectors),  and  (c,  C)  = 
vorticity  field  function  (vectors). 

Additional  notation  is  q,  =  wavenumber  variables,  qt  =  wavenumber  unit 
vectors,  and  q  =  wavenumber  vector  =  q,qv 

A  general  property  of  the  velocity  distributions  associated  with  incom¬ 
pressible  flow  is  that  the  divergence  is  zero,  as  in 

V  •  v  -  d,v,  =  0  •  (12) 

It  is  easy  to  verify  by  direct  application  of  the  divergence  operator  that  the 
divergence  of  the  vorticity  is  also  zero,  as  in 

V  c  =  did  =  0  .  (13) 

The  FT  of  equations  (12)  and  (13)  give  the  following  result: 


q,v<  =  o 

<7,0  =  0. 


(14) 


The  FT  of  the  velocity,  vorticity,  and  second-term  functions  from  equations 
(7),  (8),  and  (10)  are  as  follows: 


V0(q)  =  jn3/2[-q1q2+  q^expl-q2  /  4]/2 

C0(q )  =  n^hqM  -  q2W 3  +  %[£  +  q22]}exp[-q2 / 4]/2  (15) 

Soiq)  =  0r3/2  )1/2{-^<72+  q2q^}q3exp[-q2  /  8]/ 8 . 

In  the  FT  domain,  the  equations  giving  the  vorticity  in  terms  of  the  velocity 
are  as  follows: 

Ci  =  j{q2v 3  -  q3v2) 

Ci  =  j(q3v i  -  3)  (16) 

C3  =  j(q  1V2  ~  q2Vi)- 
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A  useful  relation  is  the  sum  of  the  above  three  components: 

Ci  +  C2  +  C3  =  j[(q 3  —  q^Vi  +  —  q^)Vi  +  (^2  —  ^1)^3]  •  (17) 

Combining  equations  (14),  (16),  and  (17)  gives  the  following  results  for  the 
V{  in  terms  of  the  C,-: 


Vi  = 


v2  = 


•f  (^2^3  ^3^2) 


(?3Ci  -  ftCs) 


,  (<7iC2  -  <?2Ci) 

3  =  T  7  ]' 

The  results  of  equation  (18)  allow  substitution  of  vorticity  components  for 
velocity  components. 

Rewriting  equation  (11)  using  the  tensor  notation  of  the  fourth  line  in 
equation  (4)  and  dropping  the  cross-product  term  results  in 

Cj(x,t)dtCj(x,t)  =  KCj(x,t)dfcj(x,t) 

d,[cj{x,  t))2  =  2KCj(x,t)dhj(x,t).  (19) 

Also  indicated  in  equation  (19)  is  the  fact  that  the  vorticity  is  a  function  of 
time.  A  useful  modification  of  equation  (19)  is  obtained  by  consideration  of 
the  second  partial  of  the  square  of  some  function,  such  as  f(x),  as  in 

d\\j{x)f  =  3x0,1 'f(x)]2}  =  3,{2 f(x)dj(x)}  =  2f(x)dlf(x)  +  2[dJ(x)f 
f(x)dlf(x)  =  (l/2)dl[f(x)f  ~ 

The  last  line  of  equation  (19)  is  rewritten  for  the  first  component  using  the 
last  line  of  equation  (20),  as  in 

3t[ci(*,f)]2  =  K{df[ci(x,t)]2  -  20,ci(x,f)]2}  •  (21) 

Similar  expressions  may  be  written  for  the  other  components  of  c(x,t).  The 
next  step  is  to  obtain  the  FT  of  the  above  expression  through  the  equation 

©o  <*>  oo  ©o  oo  oo 

3,  |  J  \{dq')3Ci(q',t)Ci(q-  q',t)  =  -Kq?  j  J  j  (dq')3ci(q',t)ci(q  -  q',t)  + 

oo  oo  oo 

2K  f  f  f  (dq'ftifflq,  -  qft]Ci{q',t)Ci(q  -  q'A) 

J  J  J  \ 


j  J  \(dq'f[dt+  Kq>-  2K(q:)(q,  -  qr)]Ci(q',t)Ci(q  -  q'J)  = 


0. 


Equation  (22)  is  solved  in  appendix  B.  The  symbol  e(x,t)  will  be  used  for 
the  enstrophy  in 

e(x,t)  =  [cj{x,t)f 
9{e(x,t)}  =  E(q,t)  =  j  J 

—  oo  —  oc 

From  appendix  B,  the  IFT  of  the  last  expression  in  equation  (23)  may  now 
be  written  as 


](dq')3Cj{q',t)Ci(q-  q',t).  (23) 


e(x,t) 


r4exp[-2x7/(4b)] 
v  (4b)7 


\{x\xz)2  +  ( *2*3 )2  +  (4b  -  xj  -  X7)2]  . 


(24) 


Equation  (24)  is  the  most  significant  and  novel  result  in  this  report.  Nota¬ 
tion  has  been  simplified  by  the  substitution  b  =  1/4  +  Kt.  As  an  elementary 
check  of  mathematical  manipulations,  it  is  easy  to  see  from  the  above  that 
the  initial  (f  =  0)  enstrophy  distribution  is  indeed  the  square  of  c0(x)  from 
equation  (8). 


5.  Confirmation  of  Results 

The  purpose  of  this  section  is  to  confirm  that  the  results  derived  in  section 
4  are,  indeed,  a  solution  of  the  transport  equation  for  instantaneous 
enstrophy  (eq  (11)).  The  original  differential  equation  (eq  (4))  is  repeated 
here: 

dtc  =  Vx(uxc)  +  KV2c  .  (25) 

Equation  (25)  is  actually  three  equations,  wherein  the  operations  are  on  the 
components  of  the  vorticity.  Again  using  the  symbol  s  for  the  second  term, 
the  three  equations  can  be  represented  as  one  equation  in  tensor  notation, 
using  the  coordinate  unit  vectors 

dtCjXj  =  Sjij  +  KQj )2CjXj  •  (26) 

The  dot  product  was  applied  to  this  vector  equation,  resulting  in  equation 
(11),  which  is  reproduced  here: 


CjdtCj  =  CjSj  +  Kcj(di)2Cj  ■  (27) 

It  is  also  necessary  to  verify  that  our  solution  satisfies  the  divergence  equa¬ 
tions  (12)  and  (13).  To  perform  this  verification,  the  solution  will  consist  of 
trial  expressions,  which  will  be  substituted  into  equation  (12),  (13),  and 
(27).  Trial  functions  will  be  designated  by  bold  letters. 

The  trial  vorticity  function  is  contained  in  the  components  of  equation  (24), 
or 


c 

L 


2exp[-x7/(4b)] 
(4b)7  7  2 
1/4  +  Kt  . 


\ 

[X1X1X3  +  X2X2X3  +  X3(4b  “  Xl  -  X2)] 

J 


(28) 
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V 

V  xv 
2A/(4L) 
v 


Note  that  the  divergence  of  the  above  is  zero  as  can  be  seen  by  summing 
the  three  components  in  equation  (29): 


3ici 


02C2 


03C3 


r2exp[-xf/(4L)] 

(4L)7/2 

^2exp[-x2/(4L)] 

(4L)7/2 

2exp[-x?/(4L)] 

(4L)7/2 


[1-  2x?/(4L)]x3 

/ 

"V 

[1—  2x^/(4L)]x3 

J 

[-2(4 L  -  x\  -  xi)x3/(4h)] 

J 


(29) 


There  is  a  direct  link  between  the  trial  vorticity  (eq  (28))  and  the  trial  veloc¬ 
ity:  find  the  FT  of  c;  apply  equation  (18)  to  find  the  FT  of  v;  and  find  the 
inverse  Fourier  transform  (IFT)  of  the  result.  This  has  not  yet  been  done. 
An  alternate  route  is  to  propose  a  functional  form  in  analogy  with  equa¬ 
tions  (7)  and  (28),  including  a  multiplicative  factor;  find  the  curl  of  that 
form;  and  derive  the  specific  form  of  the  factor  to  make  the  result  equiva¬ 
lent  to  equation  (28).  The  proposed  form  is  in  equation  (30),  where  A  is  the 
ab  initio  unknown  factor,  as  in 


A(-xix2  +  X2  xi)exp[-  xf  /  (4L)] 

[2A/(4L)][jcjxiX3  +  X2X2X3  +  X3(4F  —  xf  —  X2)]exp[~  xj  /  (4L)] 

2/(4  L)772;  A  =  1/(4L)5/2  (30) 

[1/(4L)5/2](-xiX2  +  £2 xi) exp [- xf  / (4L)]  . 

It  is  obvious  that  the  divergence  of  v  is  zero  because  of  the  opposite  sign  of 
the  two  terms,  and  the  fact  that  the  derivatives  of  the  exponential  factor 
will  make  their  magnitudes  the  same. 

To  record  the  verification  process  cogently,  set  up  the  Ac  function,  compris¬ 
ing  all  three  separately  identified  terms  of  equation  (27),  shifted  to  the  left- 
hand  side  as  in 


Ac  =  Oc  -  Pc  -  K<Pc 

Oc  —  CjdtCj 

Pc  =  C/S;  (31) 

<pc  =  C/0,)2C;  • 

Showing  Ac  is  zero  will  prove  that  the  (c,  v)  functions  from  equations  (28) 
and  (30)  are  a  consistent  solution  set  for  equation  (27).  Expansion  and  sim¬ 
plification  of  functions  (aa  (Ja  <pc)  are  done  in  appendix  C.  The  results  are 
that  pc  =  0  and  ac  =  K  <pc.  Therefore,  the  conclusion  is  that  (c,  v)  are  a  con¬ 
sistent  solution  set  for  equation  (27). 
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6.  Discussion  of  Results 


Displaying  the  nature  of  the  enstrophy  field  defined  by  equation  (24)  is 
somewhat  of  a  challenge  because  of  the  violent  variability  with  respect  to 
time.  The  size/time  scales  are  incorporated  in  the  constant  K,  so  that  any 
curves  in  terms  of  the  variables  ( x ,  t/K)  are  universal.  If  t  is  1/(4 K),  the 
enstrophy  at  the  origin  ( x  =  0)  is  already  1/32  of  the  value  for  t  =  0.  Ordi¬ 
narily,  a  time  constant  is  defined  when  some  energy  like  quantity  has  de¬ 
clined  to  one-half  of  its  starting  value.  This  would  translate  into  the  t  =  1/ 
(4 K)  if  the  numerator/ denominator  combination  was  L-1.  The  fact  that  the 
numerator/denominator  combination  changes  when  x  *  0  further  compli¬ 
cates  the  search  for  a  simple  describer  for  the  time  history.  If  the  size  is  de¬ 
fined  to  be  the  radius  of  the  circle  in  the  (xy  x2)  plane,  where  the  exponent 
of  envelope  function  is  -2,  then  the  size  a{t)  =  (4 L)1/2.  The  meaning  of  this 
is  that  the  turbule  gets  larger  with  time,  a  rather  novel  concept.  This  appar¬ 
ent  increase  in  size  is  masked  by  the  5th  power  fall  off  of  the  central  maxi¬ 
mum  of  the  distribution.  From  an  acoustic  scattering  point  of  view,  the 
scattering  cross  section  declines  with  time,  but  the  differential  scattering 
cross  section  becomes  more  peaked  in  the  forward  direction.  The 
enstrophy  is  zero  at  certain  points  in  the  {x\  —  xo)  plane  reflecting  the  fact 
that  the  vorticity  changes  sign  for  this  particular  starting  velocity  distribu¬ 
tion  as  the  point  of  interest  recedes  from  the  axis. 

To  aid  intuitive  understanding,  the  dimensionless  fields  and  variables  will 
be  discarded  and  the  pertinent  equations  will  be  as  follows: 


%,T) 


u{y,  t) 
e(y,x) 


zv(y,r) 


2  fl5  exp[-yf/(g2  +  4vt)] 


(a2  +  4vr) 


7/2 


[yi3/3y3  +  y2y2y3  + 


y3(a2  +  4vt  -  y\  -  y\)] 
Q0fl5exp[-t/f/(fl2  +  4vr)] 


( a 2  +  4vt  )' 


5/2 


[—  3/ 1  y  2 +  y2y^ 


^4fl2g10exp[-2i/- /(g2  +  4vt)] 
v  ( a 2  +  4vt)7  j 

2  ..2x2 


[(y3y3)2  +  (y2y3)2  + 


(a2  +  4vt  -  y1  -  y2)  ] 

pQofl10exp[-2i/f/(fl2  +  4vt)] 
2  (a2  +  4vt)5 


[y2  +  y\]- 


(32) 


Perhaps  the  best  way  to  gain  some  numerical  idea  as  to  the  variations  of 
the  turbule  enstrophy  is  to  consider  the  total  of  this  quantity.  As  a  matter  of 
convenience  this  is  done  in  the  following  equation,  where  the  volume  inte¬ 
gral  of  the  enstrophy  is  calculated  as: 


*(T)  =  J  j  j  ^yfe{y,x) 


5,r3/W 0  ) 

4^2 (a2  +  4vt)7/2J  ‘ 


(33) 


9 


The  new  function  w(y,  t)  in  equation  (32)  is  the  kinetic  energy  concentra¬ 
tion.  The  total  kinetic  energy  is  calculated  from  the  following  expression: 


GX5  ©O  ©O 

W(t)  =  j  J  J  {dyfw{y,  t) 


W(t) 

W(  0) 


(4L)_5/2 . 


jt3/W°P 
842(a2  +  4vt)5/2 


(34) 


From  the  last  expression  in  equation  (34),  let  f1/2  be  that  time  for  which  the 
ratio  is  1/2  (a  condition  that  occurs  when  Kt\/2  =  (22^5  - 1)/4  =  0.0798770). 
The  intent  is  now  to  investigate  the  state  for  seven  different  turbule  sizes. 
The  relevant  data  is  collected  in  table  1 . 


The  sizes  are  from  the  ensemble  reported  previously  [1].  In  this  ensemble, 
the  characteristic  size  of  the  largest-sized  turbule  was  a\  =  10  m.  The  maxi¬ 
mum  velocity  in  the  largest  turbule  was  approximately  1  / 100  of  the  speed 
of  sound  or  V\  =  3.44  m  •  s-1.  The  turbule  size  ratio  was  an/an+ \  =  0.827787. 
For  a  velocity  structure  constant  of  0.1111  m4//3  •  s-2,  the  turbule  spacing 
parameter  was  €  =  8.  From  equation  (6),  the  maximum  velocity  occurs  for 
y3  =  0  on  a  circle  with  a  radius  of  pm  =  2~1/2  a.  The  kinematic  viscosity  coef¬ 
ficient  is  taken  to  be  0.15  x  10-4  m2  •  s"1.  The  starting  maximum  velocity  of 
the  different-sized  turbules  is  assumed  to  vary  according  to  the  (1/3) 
power  rule  applied  to  the  size  ratio.  The  critical  Reynolds  number  for 
boundary  layer  flow  with  zero  pressure  gradient  is  given  [7]  as  about  600, 
indicating  the  value  when  the  flow  becomes  unstable.  The  sizes  in  the  table 
cover  a  Reynolds  number  range  spanning  this  value. 

Other  valuable  information  is  provided  in  table  2,  including  data  on  total 
kinetic  energy  and  the  pressure  differential  A p  at  the  origin,  which  is  calcu¬ 
lated  from  the  following  relationship  to  the  velocity  in  the  y\  -  y2  plane: 


rdp 

Ap(r) 


p(u2  /  r)rdr  r 

2  J  drzv(r,  0,t)/ r 
o 


>/y?  +  y\ 
r  qW^p  ' 

^4  (a2  +  4vt)4, 


(35) 


Table  1.  Turbule  decay  data. 


Characteristic 
size,  a 
(m) 

Maximum 
velocity,  vm 
(m  •  s"1) 

Radius  for 
maximum,  pm 
(m) 

Angular 

velocity, 

Qq 

(s-1) 

Decay 

constant, 

K 

Relative 
half¬ 
time,  t\/2 

Half¬ 
time,  Tj/2 
(s) 

Enstrophy  int., 
«(y/Ti/2) 
(m3  •  s-2) 

1.04  x  10 

1.62x10 

7.32  xlO-1 

1.63  x  101 

8.58  x  lO"7 

9.30  x  104 

5.71  x  103 

3.10  xlO2 

3.33  x  10-1 

1.11x10 

2.36  x  10"1 

3.47  xlO1 

3.89  x  10-6 

2.05  x  104 

5.91  x  10 2 

6.19  x  101 

1.07  x  10-1 

7.58  x  10”1 

7.58  x  lO”2 

7.40  x  101 

1.77  xlO"5 

4.52  x  103 

6.12  xlO1 

9.35  x  10 

3.45  x  10“2 

5.20  x  lO”1 

2.44  x  10“2 

1.58  xlO2 

8.01  x  10"5 

9.97  x  102 

6.33  x  10 

1.41  x  10 

1.11  x  lO"2 

3.56  x  lO"1 

7.84  x  lO-3 

3.35  x  102 

3.63  x  10-4 

2.20  x  102 

6.55  x  lO"1 

1.79  x  10”1 

3.57  x  10-3 

2.44  x  10-1 

2.52x10-3 

7.14  x  102 

1.65x10-3 

4.85  x  101 

6.78  x  10-2 

1.79  xlO'1 

1.15  x  10-3 

1.67  xlO-1 

8.12  xlO-4 

1.52  xlO3 

7.48  x  lO-3 

1.07  xlO1 

7.02  x  lO"3 

1.79  x  lO"1 

Note:  v=2.50  x  10~ 5  m2  *  s  J  and  K  tyi  -  0,079877. 
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The  maximum  pressure  differential  (in  the  largest  turbule)  is  about  86  Pa, 
compared  to  the  nominal  atmospheric  pressure  of  101,000  Pa.  These  pres¬ 
sure  data  confirm  that  the  constant  density  assumption  is,  essentially, 
valid. 

A  comparison  between  the  decay  time  of  this  theory  and  the  decay  time 
predicted  by  Kolmogorov's  theory  is  instructive.  The  relation  from  which 
the  latter  time  may  be  calculated  is  as  follows: 

u  -  (a2/e)'/3  ■  (36) 

In  equation  (36),  £  is  the  energy-decay  rate  per  unit  mass.  A  typical  value 
[8]  for  it  is  £  =  10'3  m2  ■  s'3.  Table  3  contains  the  relative  comparative  data. 

From  table  3,  we  see  that  the  comparable  times  occur  when  the  size  is 
around  1  cm.  The  theory  of  this  report  shows  a  steeper  decline  for  smaller 
sizes  than  Kolmogorov's  prediction,  suggesting  that  the  present  theory 
might  be  useful  for  analyzing  the  turbule  interaction  in  this  region.  For 
example,  consider  two  turbules  with  initial  velocity  distributions  similar  in 
form  to  that  considered  in  the  table,  separated  by  a  suitable  distance.  The 
only  terms  in  equation  (11)  initially  present  would  be  the  interaction  terms, 
and  the  decay  (or  growth)  of  the  interaction  could  prove  to  be  very 
interesting. 


Table  2.  Ancillary  turbule  data. 


Characteristic 

Maximum 

Radius  for 

Angular 

Kinetic 

Half-kin. 

Center 

Half-center 

size,  a 

velocity,  vm 

maximum,  pm 

velocity. 

energy. 

energy 

pressure 

pressure 

(m) 

(m  •  s_1) 

(m) 

Qq 

w 

Wj/2 

Ap 

Api/2 

(s'1) 

a) 

0) 

(Pa) 

(Pa) 

1.04x10 

1.62x10 

7.32  x  10-1 

1.63  x  101 

1.88  xlO2 

9.41  x  101 

8.62  x  101 

2.84  x  101 

3.33  xlO"1 

1.11x10 

2.36  xlO'1 

3.47  x  101 

2.94  x  10 

1.79x10 

4.05  x  101 

1.33  x  101 

1.07  xlO"1 

7.58  x  10-1 

7.58  x  lO"2 

7.40  x  101 

4.60  x  lO'2 

2.81  x  lO’2 

1.90  xlO1 

6.27  x  10 

3.45  x  10-2 

5.20  x  lO"1 

2.44  x  lO"2 

1.58  x  102 

7.20  x  10^ 

4.39  x  lO"4 

8.92  x  10 

2.94  x  10 

1.11  x  10-2 

3.56  x  10"1 

7.84  x  lO"3 

3.35  x  102 

1.13  xlO'5 

6.86  x  lO"6 

4.19  x  10 

1.38  x  10 

3.57  xlO"3 

2.44  x  10"1 

2.52  x  lO'3 

7.14  x  102 

1.76  xlO'7 

1.07  xlO"7 

1.97x10 

6.49  x  lO'1 

1.15  xlO”3 

1.67  xlO"1 

8.12  xlO"4 

1.52  xlO3 

2.75  x  lO'9 

1.68  xlO"9 

9.24  x  lO"1 

3.05  x  10"1 

p0=1.01  xl05Pa  p=  1.21  kg  -m  3 


Table  3.  Decay  time  comparison  with  Kolmogorov  prediction. 

a  (m)  1.04x10  3.33  xlO"1  1.07  xlO'1  3.45  x  10'2  1.11  x  HT2  3.57  x  KT3  1.15  xlO'3 

xK(s)  1.03  xlO1  4.80x10  2.25x10  1.06x10  4.98  x  10'1  2.34  xlO"1  1.10  xlO'1 

x1/2(s)  5.71  xlO3  5.91  xlO2  6.12  xlO1  6.33x10  6.55  xlO'1  6.78  x  IQ-2  7.02  xlO'3 
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7.  Conclusion 

The  mathematical  developments  reported  in  this  document  contain  a  valid 
solution  to  the  enstrophy  correlate  of  the  N-S  equation,  giving  the  time  his¬ 
tory  of  the  decay  of  the  enstrophy  from  a  particular  choice  of  the  initial 
velocity  distribution  of  an  isolated  velocity  turbule.  There  may  be  other 
initial  velocity  distributions  that  lead  to  a  similar  analytical  solution. 

We  found  that  the  time-dependent  velocity  field  derived  here  is  not  a  solu¬ 
tion  to  equation  (2),  the  N-S  equation.  Details  of  this  analysis  are  given  in 
appendix  D. 
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Appendix  A. — Applicable  Fourier  Transform  Theory 

Fourier  transform  (FT)  theory  was  used  in  section  4  of  the  main  report  to 
solve  the  equation  resulting  from  forming  the  dot  product  of  the  vorticity 
and  equation  (3).  The  applicable  theory  is  outlined  in  this  appendix,  along 
with  the  notation  employed. 

The  FT  of  function  g(y)  of  variable  y  to  variable  /is  defined  as  follows: 


=  G(f)  =  J  dyg(y)exp(jfy)  .  (A-l) 

— oo 

The  inverse  FT  (IFT)  of  function  Gif')  of  variable  f'  to  variable  y '  is  defined 
as  follows: 

®-l{G(f'),y'}  =  g(y')  =  (2k)-1]  df'  G(f')  exp(-;/'y')  .  (A-2) 

—  OO 

The  Dirac  delta  function,  the  symbol  for  which  is  8(y),  and  its  FT  are  of 
considerable  importance.  8(y)  is  defined  by  the  limiting  process  recorded 
below: 


S(y)  =  lim  —I;  -A/2  <  y  <A/ 2-0 

A  ->  0  V-d. 


=  0;  y  <-A/2  ox  A/2<y- 

The  integral  of  8 (y)  is,  therefore,  unity.  The  FT  of  8(y)  is  the  following: 

°°  /  i  \ 

9{5(y))  =  f  dy  8(y)  exp(;/y)  =  lim  —  f  dy  exp(jfy)  = 

J-  A-toiAj^  2 


(A-3) 


lim 

A  ->  0 


exp(//y) 

A/2 

=  lim 

2sin(d//2 )' 

;a/  J 

-A//2 

A  — >  0 

L  Af  J 

=  1  • 


(A-4) 


The  FT  of  8(y)  is,  therefore,  a  constant;  i.e.,  it  has  a  uniform  spectrum.  The 
IFT  of  1,  then,  recovers  the  delta  function 


yMl)  =  (2k)-1  J  df  exp (-)fy')  =  S(y')  _  (A.5) 

— oo 

When  this  integral  construction  is  encountered,  substitution  of  a  suitable 
delta  function  is  appropriate. 

The  next  few  expressions  define  the  convolution  integral,  which  is  useful 
when  dealing  with  the  product  of  two  functions. 

f(y)  =  g(y)Ky) 

oo 

=  F(f)  =  J  dyg(y)h(y)exp(jfy)  •  (A-6) 

—  oo 
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In  equation  (A-6),/(y)  is  defined  as  the  product  of  the  functions  giy)  and 
h(y).  The  FT  of/(y)  is  defined  the  standard  way.  To  develop  the  convolution 
integral,  giy)  and  h(y)  are  expressed  in  terms  of  the  IFT  of  their  transforms 
Gif)  and  Hif),  as  in 

oo  ©o  00 

J  dy  (2tt)_1  J  df'G(f')exp(-jf'y)(2K)~l  J  df'H(f")expi-jf"y)expijfy) 

OO  00  ©o 

(2 ny2  J  df'Gif')  J  df'Hif")  J  dy  exp[(j  y)(f  -  f  -  /")] 

00  ®o 

(2a:)-1  J  df'Gif)  J  df'Hif")8if-  f  -  f")  (A-7) 

«—  ©C  ^OO 

OO 

(2/r)-1  \df'Gif)Hif-  f)  =  (2 nf  Gif)* Hif)  . 

~OG 

The  last  line  of  equation  (A-7)  defines  the  convolution  integral.  The  right¬ 
most  member  of  the  last  line  shows  how  the  integral  is  sometimes  repre¬ 
sented.  The  convolution  integral  is  unchanged  by  interchange  of  the  func¬ 
tional  dependencies  between  G  and  H. 

The  three-dimensional  FT  of  function  hix)  of  vector  variable  x  to  vector 
variable  q  is  defined  as 


OOOOOO 


9{hix),q}  =  J  J  J  idx)3hix)expijqjXi)  =  Hiq)  (A.g) 


—OO  — 00  —00 


The  IFT  of  function  Hiq')  of  vector  variable  q  to  vector  variable  x'  is  shown 
as 


l{Hiq'),x'\  =  i2n)-3  J  J  \  (dx'f  Hiq')  exp(-j q'X.)  =  hix')  . 


(A-9) 


The  FT  of  the  partial  derivative  of  a  function  (provided  giy)  approaches 
zero  sufficiently  rapidly  at  large  y)  is 


jr g(y)'f 
dy 


=  -if  Gif) 


Other  notation  conventions  were  included  in  the  text  of  this  report. 


(A-10) 


Appendix  B. — Equation  Solution  Procedure 

The  Fourier  transform  (FT)  shown  in  equation  (22)  of  this  report,  is  repro¬ 
duced  here  in  the  equation  below. 


J 


1 


oo  oo  oo 


W)3Ci(q',t)Ci{q-  q'.t)  =  -Kq]  j  J  J  (dq')3Ci(q',t)Ci(q  -  q’A)  + 

— -  oo  — oo  — oo 

oo  oo  oo 

2Kj  J  tf.t) 


oo  oo 


j  J  j(dq'f[dt+  Kqf-  2K(q')(qi-  q^]Ci(q',t)Ci(q  -  q'J)  =  0. 


(B-l) 


The  solution  to  equation  (B-l)  is  carried  out  in  equation  (B-2),  as  well  as 
conversion  back  into  the  space  domain. 

To  ensure  that  the  integral  in  the  second  line  of  equation  (B-l)  is  zero,  let 
the  integrand  be  zero.  Temporarily  indicating  the  vorticity  product  on  the 
right  as  the  function  P(q,  q',  t),  the  Laplace  transform  (LT)  of  the  integrand 
may  be  calculated  as 


#{P(<M',f)}  =  tfKq,q\s)  =  jdtP{q,q'rt)exp(-st) 


[a,  +  K(q,)2  -  2K{q’i){qi  -  q;)]P(q,q',t)  =  0 
[s  +  K( q;)2  -  2 K(q'i)(qi  -  q;)Mq,q',s)  =  P(q,q',0) 

P{W’S)  [S+  Riqf-  2K{q'i){qi  -  <?/)] 
P{q,q',t)  =  P(q,q'r0)exp{-K[(qj)2  -  2(q'i)(qj  -  ql)]t} . 


(B-2) 


The  first  line  in  equation  B-2  defines  the  LT.  The  second  line  is  the  inte¬ 
grand  from  equation  (B-l).  The  third  and  fifth  lines  of  the  equation  are 
standard  LT  [6].  The  left-hand  side  of  the  last  line  of  the  equation,  minus 
the  right-hand  side  may  now  be  substituted  for  the  integrand  of  the  last 
line  of  equation  (B-l),  as  in 


J  J  J  ( dq'?Ci(q',t)Ci(q  ~  q'.t)  = 

— oo  — oo  — oo 

oo  oo  oo 

j  J  j(dq'fUq',0)Uq-  q'r0)exp{-K[(qt)2-  2(q;)(qi-  q;))t}. 
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gj-MLHS) 


The  next  step  involves  the  IFT  of  equation  (23).  Addressing  the  left-hand 
side,  the  IFT  integration  will  be  applied,  and  the  integral  expression  for  the 
C's  will  be  substituted,  as  in 


(2 tt)'3  j  j  “j  (dq)3exp(rjqixd’j  j  j  (dq')3Ci(q',t)C\(q  -  q',t) 

€>©  OO  ©O  ©O  ©O  ©O 

(27T)"3  J  J  $  (dq)3exp(-jqiXi)  J  J  J  (dq'f 

_  oo  —  OO  —  oo  -®o  — oo  -oo 

j  j  ~j{dx')3exp(jqlxl)ci(x'tt) 
j  j  j  (dx")3exp[/0?,  -  ql)x']ci(x",t) 

—  OO  —  ©O  — ©Q 

©o  ©o  ©o  ©o  ©o  ©o 

(2tt)“3  J  J  ](dq')3\  J  \  (dx')3exp(jq'iX;)ci(x',t) 

—  OO  —  OO  — OO  — OO  — OO  — oo 

©O  GO  ©O 

J  J  J  (dx"fexp[-jq'iX']ci(x",t) 

J  }  J  (dq)3exp[jqi(x''-X,)] 

—  OO  —  OO  —oo 

oo  ©o  oo  ©o  oo  ©o 

J  J  J  ( dq'f  J  J  J  (dx'fexpijq  ,'Xj) Ci (x'r t) exp[-jq'j  Xi]ci(x,  t ) 

j  J  ]  (dX'fci(x',t)ci(x,t)  j  j  j  (dq'fexp[Jq;(x;  -  *,•)] 

(2^r3  [ci(x,o]2. 

The  square  of  the  vorticity  was  defined  to  be  the  enstrophy  in  section  4 

e(x,t)  =  [cj(x,t)f 

&{e(x,t)}  =  E(q,t )  =  (27r)~3  J  J  J  ( dq')3Cj(q',t)Cj(q  -  q'J) 


(B-4) 


E(q,t)  =  (2 TT)"3  j  J  j(dq')3P(q,q',t). 


(B-5) 
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oo  oo  oo 

j  I  Jw3  (?!“  <)(^3-  ^)expKfl,  -  <?,')2/4] 

— oo  — oo  —  oo 

exp [~K(qi  -  q'i)2t]exp[-j(qi  -  q')Xi] 

oo  oo  oc 

e2(ac,0  =  (tt3/4)  J  J  j  (dq')3q2q'3exp[-{q')2/A]exp[-K(q;)2t]exp(-jq'iXi) 


J  j  j  2  -  <fe)(93  -  '?3)exp[-(<?,  -  qi)2/  4] 

—  OO  —  oo  — oo 

exp[-K(q,  -  q')2t]exp[-j(q,  -  q',)x,} 
e3(x,t)  =  (n3/ 4)J  j  ]  (dq'fKqi)2  +  (^)2]exp[-(^)2/4] 

—  oo  — oo  — oo 

exp[-K(g,'  )2  f]exp(-;<7,'  *,) 

J  j  J(^)3[(^i-  <)2+(^2-  <?2  )2]exp[-(<7,  -  <?/ )2/4] 

— oo  —  oo  — oo 

exp[-K(<j,.  -  q'i)2t]exp[-j(qj-  q',)x,]- 

The  right-most  integrals  are  IFTs  and  produce  functions  that  have  no  q 
dependence.  The  notation  has  been  simplified  by  setting  L  =  1/4  +  Kt. 
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ei(x,t)  =  (tt3/4)  J  j  J(^')3^'exPH<7/)2/4]exp[-K(^)2i]exp(-/(?;x,) 


(-xlX3exp[-xim)]\ 

l  32  L7/V/2  J 

e2(*,0  =  (tt3/4)  J  j  j  (^,)3^2'?3exPK^')2/4]exp[-K(^')2i]exp(-^'x1) 


^-X2X3exp[-^/(4L)]^ 

l  32  L7/V/2  J 

e3(^f)  =  (tt3/4)J  j  ]  (dq'f[(cj2)2  +  (^)2]exp[-(^')2/4] 


(B-8) 


exp[-X(^'  )2  f]  expo's, ) 


^(4 L-  x2  -  x2)exp[-x?/(4L)p 
,  32  L7/V/2  ; 


The  integrals  in  equation  (B-8)  are  the  same  as  the  right-most  integrals  of 
equation  (B-7),  so  that  the  result  simply  squares  the  (  )  brackets. 


e{x,t)  = 


^4exp[-2x2/(4L)] 

v  m7 


[(xiX3)2  +  (X2X3)2  +  (4T  -  x\  -  X2)2] 

J 


(B-9) 


Equation  (B-9)  is  the  formal  result  of  this  investigation,  and  was  reported 
in  section  4  as  equation  (24). 
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Appendix  C. — Details  of  the  Confirmation  Process 

The  purpose  of  this  appendix  is  to  record  the  details  confirming  that  the 
enstrophy  expression,  equation  (24)  of  the  main  report,  is  a  solution  of  the 
transport  equation  for  instantaneous  enstrophy,  equation  (27).  In  section  5, 
trial  vorticity  and  trial  velocity  functions  (denoted  by  bold  letters  c  and  v, 
respectively)  were  inferred  from  equation  (24)  and  are  reproduced  here  in 
equations  (C-l)  and  (C-2). 


c 


2exp[-x2/(4  L)] 
(4L)7/2 


\ 

[*1*1X3  +  *2*2*3  +  *3(4T  -*i  -*2)]  . 

J 


(C-l) 


v  =  [1/(4L)5/2](-x1x2  +  *2 *1) exp[- xj / (4L)]  •  (C-2) 

I  showed  in  section  5  that  the  divergence  of  c  and  v  are  zero — a  condition 
established  early  in  this  investigation.  The  confirmation  process  that  began 
in  section  5  identified  a  function  Ac  made  up  of  the  three  terms  of  equation 
(27).  These  three  terms  were  separately  identified  in  equation  (31),  which  is 
reproduced  here  as  equation  (C-3). 

Ac  =  etc  ~  Pc  -  *<Pc 
ode  =  CjBtCj 

Pc  =  CjSj  (C-3) 

<pc  =  Cj(di)2Cj  ■ 

Expansion  of  function  ac  is  formulated  in  equation  (C-4). 


Ode  =  Cl3tCl  +  C23tC2  +  C33tC3 

=  2  xi  *3  exp[-  x2  /(4L)](4L  )~7 /20,{2  xi  *3exp[-  xj  /  (4L)](4L  )~7/2} 

+  2x2X3exp[-xf/mmr7/2dl{2x2x3exp{-xf/mmr7/2} 

+  2(4 L-  x}~  *2)exp[-x?/(4L)](4L)”7/2 
3,{2(4L  -  xl-  *2)exp[-x?/(4L)](4L)_7/2} . 

Taking  the  time  derivatives  and  simplifying  gives  the  following: 

ac  =  [2  xi  *3  )2exp[—  xf  /  (4L)](4L  )~7 1 2  {(4K)  exp[-  xj  /  (4L)](4L  )'u  /  2 
[xf-  (7/2)(4L)]}  +[2x2*3]2exp[-x?/(4L)](4L)-7/2{(4K) 
exp[-x2/(4L)](4L)-n/2[x2-  (7/2)(4L)]}  +  4(4L  -  x?  -  xl) 
exp[-x2/(4L)](4L)-7/2{(4K)exp[-x?/(4L)](4L)-11/2 
[(-7/2)(4L  -  xl  -  *i)(4 L)  +  x2(4L  -  xf  -  xl)  +  (4L)2]}  . 
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Vc 


Even  further  simplification  yields  the  following: 

etc  =  (8fC)exp[-2x2/(4L)](4L)'9{x2(xi  +  x2)[2x2  -  7(4L)] 

+  2x2(4L-  x\  -  xlf  -  7(4L)(4L  -  x\  -  x2)2  (C- 6) 

+  2(4L)2(4L  -  x\  -  *?)}• 

Expansion  of  function  <pc  is  shown  in  equation  (C-7): 

Cl5?Cl  +  C2djci  +  C30?C3 

2  xi  xj  exp[-  x?  /  (4L)](4L  )-7/  2  a2  {2  xi  xz  exp[-  xf  /  (4L)](4L  )~7/2\ 

+  2x2X3exp[-x?/(4L)](4L)-7/23?{2x2X3exp[-x?/(4L)](4L)-7/2} 

+  2(4L  -  xi  -  xi)exp[-x2/(4L)](4L)-7/2  (C'7) 

x  32{2(4L  -  xi  -  ^)exp[-x?/(4L)](4L)-7/2}. 

Breaking  equation  (C-7)  into  parts  and  performing  the  differentiations 
gives  the  following  three  equations. 


<Pcl  =  4xix2exp[-x2/(4L)](4L)'73?{xiexp[-x?/(4L)]) 

+  4*2*3  exp[-  xf  /  (4L)](4L  )‘7  di  {exp[-  x2  /(4L)]) 

+  4(4L  -  xi  -  *2)  exp[-  xf  / (4L)](4L  )-7  a2 
{(4L  -  xi  -  x2)exp[-x2/(4L)]) 

=  8x2xi[2x2 /(4L)  -  3]exp[-2x?/(4L)](4L)-8 
+  8x2X3 [2 xi /(4L)  -  l]exp[-2x2/(4L)](4L)-8 
+  8(4L-  xi  -  xlWL-  xi  -  xi)(2x?/(4L)  -  1)  +  2xfl 
exp[-2x?/(4L)](4L)-8. 


<Pc2  =  4[xix3]2exp[-x2/(4L)](4L)'7al{exp[-x2/(4L)]) 

+  4x2x1  exp[-  xf  /(4L)](4L  )'7  a2  {x2  exp[-  xf  /  (4L)]} 

+  4(4L  -  xi  -  X2)  exp[-  xf  / (4L)](4L  )~7 dj 
{(4 L-  xi  -  x2)exp[-x2/(4L)]} 

=  8[xix3]2[2x2/(4L)  -  l]exp[-2x2/(4L)](4L)-8 
+  8[x2x3]2[2x2/(4L)  -  3]exp[-2x2/(4L)](4L)"® 

+  8(4L-  xi  -  X2)[(8E  -  xi  -  x2)(2x2/(4L)  -  1)  +  2x2] 
exp[-2x2/(4L)](4L)-8. 
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<pc3  =  4xh3exP[-xUmmr7di{x3exP[-xf/(m) 

+  4x1  x3exp[-  x} /(4L)](4L  r7  a  1  (x3exp[-  xl  /  (4L)]} 

+  4(4 L-  x2  -  xl)exp[-xl/(4L)](4L)“733 
{(4L  -  x\  -  xl)exp[-xl/(4L)]} 

=  8x?x§[2xl/(4L)  -  3]exp[-2^?/(4L)](4L)-8  (C-10) 

+  8xlxl[2x23/(4L)  -  3]exp[-2x?/(4L)](4L)-8 
+  8[2x3/(4L)  -  1](4L  -  xi-  xl  ? exp[-2  xl / (4L)](4L )"* . 

Summing  up  the  expressions  in  equations  (C-8),  (C-9),  and  (C-10)  and  sim¬ 
plifying  the  results  achieves 

<PC  =  8xlxl[2xl/(4L)  -  3]exp[-2x2/(4L)](4L)-8 
+  8x2x2 [2  x2/(4L)  -  l]exp[-2x? / (4L)](4L)"® 

+  8(4L  -  xi  -  xi)[(8L  -  xi  -  X2)(2x2/(4L)  -  1)  +  2x2] 
exp[-2x2/(4L)](4L)-8  +  8[xix3]2[2xl/(4L)  -  1] 
exp[-2x?/(4L)](4L)-8  +8[x2x3]2[2xl/(4L)  -  3] 
exp[-2x?/(4L)](4L)_s  +  8(4L  -  x2  -  xl)[(8L  -  xi  -  xl) 

(2 xl  / (4L)  -  1)  +  2  x2]  exp[-2 xi / (4L)](4L  )"*  (C-ll) 

+  8x2x2[2x3/(4L)  -  3]exp[-2x?/(4L)](4L)-8 
+  8x2x3[2x3/(4L)  -  3]exp[-2x?/(4L)](4L)'8 
+  8[2x3/(4L)  -  1](4L  -  xi2  -  xl )2  exp[-2  xl  / (4L)](4L  )-8 . 

Further  simplification  yields  the  following  equations. 

<pc  =  8 exp [-2  xi  / (4L)](4L )“9{[x3 xi  f  [2  xi  -  3(4 L)]  +  [x2x3]2[2x(  -  (4L)] 

+  (4L  -  xi  -  xl)[(8 L  -  xi  -  xl)(2xi  -  (4L))  +  2(4L)x?] 

+  [xix3]2[2x2  -  (4L)]  +  [x2x3]2[2x2  -  3(4L)] 

+  (4L  -  xi  -  xl)[(8L  -  xi  -  xl)(2xl  -  (4L))  +  2(4L)x2]  (C-12) 

+  [xix3]2[2x3  -  3(4L)]  +  [x2X3]2[2x3  -  3(4L)] 

+  [2xl-  4L](4L  -  xi  -  xl)2}- 

<pc  =  8 exp[-2 xi / (4L)](4L )-9 |x3 [x2  +  xl][2xl  -  7(4L)]  +  2x1(4 L  -  xi  -  xl? 

-7(4L)(4L  -  xi  -  xl?  +  2(4L)2(4L  -  xi  -  xl))  •  “  ^ 
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Comparing  equation  (C-6)  with  equation  (C-13),  we  see  that  <*c  -  K  (pc,  so 
they  will  cancel  each  other  out  when  they  are  entered  into  the  Ac 
expression. 

The  second  term  from  equation  (C-3)  is  more  conveniently  written  in  vec¬ 
tor  notation  as 

Pc  =  c[Vx(vxc)].  (C-14) 

First  calculate  the  cross  product  of  the  velocity  and  the  vorticity: 

vxc  =  -exp[-2x2/(4L)][(4L)“V2]£3x3[xi  +x?> ]  ■  (C-15) 

Next,  find  the  curl  of  the  expression  in  equation  (C-15),  which  will  be  the 
function  s,  as  in 

s  =  exp[-2x^/(4L)][(4L)_6]x3[4b  -  x\  -  x^\[~x\X2  +  £2*1]  •  (C-16) 

The  function  Pc  will  then  be  the  dot  product  of  c  with  the  above,  or 


2exp[-3x?/(4L)] 


*3  [4 L  -  x\ 


(4L)i9/2  r  1 

[£1X1X3  +  £2X2X3  +  £3(4L-  xi  -  x2)]  [-£1X2  +  £2X1] 


Therefore,  the  second  term  is  identically  zero. 


(C-17) 
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Appendix  D. — Analysis  of  the  Enstrophy  Equation 
Applied  to  the  Navier-Stokes  Equation 

In  this  appendix,  the  inferred  time-dependent  velocity  field  will  be  applied 
to  the  Navier-Stokes  (N-S)  equation  to  determine  if  it  is  a  solution.  The 
form  of  N-S  to  be  used  is  equation  (2),  which  is  reproduced  here  as  equa¬ 
tion  (D-l). 


du 

dr 


\r  / 


I Vp  -  (1  / 2)V u2  +  MX(Vxw) 
V(V  •  u)  . 


V2u  + 

fi  + 

IpJ 

(f 

3  Py 

(D-l) 


The  trial  vorticity  and  velocity  field  equations  are  reproduced  here  from 
equation  (32). 


<%,t)  = 


u(y,T)  = 


2fl0g5exp[-yf/V  +  4vt)] 
(a2  +  4vt)7/2 

y3(fl2  +  4vt  -  y]  -  y\)\ 

Q0fl5exp[-y- /(a2  +  4vt)] 


IViViVs  +  ViViVz  + 


(a2  +  4vt) 


,5/2 


(D-2) 


l  y \V 2  +  y2yJ  ■ 


An  expression  for  the  pressure  will  be  required;  that  expression  will  be 
developed  after  the  form  of  equation  (34)  and  shall  be  shown  as  equation 
(D-3). 


r'dp  =  p(u2  /  r')r'dr'  r'  =  -Jy'l+y'l 


p(°o,r)  -  p(y,r)  =  2  J  dr'w(r'/y3,t) / r' 


&o  awp 


4(g2  +  4vt)5 

p(y,  ^)  =  pK  t)  - 


^i2+y2 

exp[-2  y2  /  (g2  +  4vt)] 


(D-3) 


4(a2  +  4vt)5 


exp[-2y-  /(a2  +  4vt)] 


To  carry  out  the  confirmation  process,  trial  expressions  will  be  identified 
by  bold  symbols.  Au  will  represent  equation  (D-l)  with  all  the  terms 
moved  to  the  left-hand  side  (except  for  the  last  term,  which  is  zero  because 
the  divergence  of  the  velocity  is  zero).  If  Au  is  zero  when  the  velocity  ex¬ 
pression  of  equation  (D-2)  is  substituted,  equation  (D-2)  is  a  solution  to 
equation  (D-l).  Two  trial  expressions  will  also  be  calculated:  Ai  and  A2. 
These  functions  represent  the  differently  orientationed  components  of 
equation  (D-l).  If  Ai  and  A2  are  simultaneously  zero,  Au  will  be  zero.  These 
identifications  are  made  in  the  equation  (D-4). 
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Au  =  Aj+A2  A^aj+Pj  +  tpj  A2  =  a2+P2 

)vP  = 


'l' 


VP) 


/  A 

tlpA10 

(a2  +4  vt)6 


P,  =  (1/2)Vu2  =  j-  D°fl 


2-10 


',y,]exp[-2y,2/(fl2+4vT)] 

v{[y2  +yi2]exp[-2yf  /(«2 +4vr)]} 


ip1  =  u  x  d  =  — 


p(a2  +  4vt)  J 

exp[-2y2/(fl2  +4vt)] 


,  2  _10 


(a2  +4vt) 

{[-  y^y  i + y2yi] x  [y^iys  +  y2y2y3  +y3(«2  +  4vT  -  y?  -  y2)] } 


(D-4) 


ol2  =  dzU  =  dx< 

p  ■> 

Q0a5 

/  0  \5 / 2 

>{[-  yiy2+y2yi]exp 

-y2/(A2+4vr)]} 

(a2  +  4vt) 

l 

P2  =  v  V2u  =  -v 


[a2  +  4vt)5/2 


V2{[-yay2  +y2yi]exp[-y,2  /(a2  +4vt)]}. 


The  gradient  operation  has  already  been  accomplished  for  04,  and  no  fur¬ 
ther  simplification  is  needed.  Next,  simplify  the  Pi  function: 


Pi  = 


<  ^a10  ^ 


{(a2  +  4vr)[y1y1  +  y2y2]  -  2[y2  +  y2] 


(a2  +  4vr  )6j 

[y2yi  +  y2ya  +  h  y3llexPl-2yf/(fl2  +  4vT ')!  • 

The  simplify  the  function: 
f  2  £io  A10  A 


(P-5) 


<Pi  = 


exp[-2y2/(a2  +  4vr)][y1y1(a2  +  4vr-  y2  -  y2) 


(D-6) 


Ua2  +  4vr)6, 

+  y2y2(«2  +  4vt  -  y2  -  y2)  -  y3y3(y? +  y2)l  • 

Find  the  sum  of  ai,  pi,  and  <pp 
Ai  =  ai  +  Pi  +  <Pi 

4vt)]\  _  (a2  +  4vr)(.  +  - „  ,  (D-7) 

(a2  +  4vt  )  J 

In  doing  so,  we  find  that  Ai  is  not  zero.  Proceed  by  simplifying  the  A2  func¬ 
tion  by  first  expanding  ot2: 


a2  = 


f  Ho  A5 


I-yjy2  +  y2yJexp[-y2/(A2  +  4vr)] 


Ua2  +  4vr)9/2y 
(4v)[y2  -  (5/2)(a2  +  4vr)]  . 


(D-8) 
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Appendix  D 


Finally,  expand  p  2: 


-4v 


Do  a 


( a 2  +  4vt  y 


[-Pi  y2  +  y2y^M-y2i^a2  +  4vr)] 


[yf  -  2(a2  +  4vt)]  . 


(D-9) 


Comparing  equation  (D-8)  with  equation  (D-9),  we  see  that  A2  is  not  zero, 
although  the  difference  has  the  time  dependence  of  the  vorticity,  but  is  in 
the  opposite  direction  of  the  velocity. 

Since  Aj  and  A2  are  each  not  zero  and  have  different  dependence  in  the 
exponential  factor,  their  sums  A  cannot  be  zero,  and  hence  equation  (D-2) 
is  not  a  solution  to  equation  (D-l). 
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